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Abstract 

The strong coupling constants of heavy baryons with light vector mesons are 
calculated in the framework of the light cone QCD sum rules using the most general 
form of the interpolating currents for the heavy baryons. It is shown that the sextet- 
sextet, sextet-antitriplet and antitriplet-antitriplet transitions are described by one 
invariant function for each class of transitions. The values of the electric and magnetic 
coupling constants for these transitions are obtained. 
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1 Introduction 



Experimental and theoretical studies of charmed and bottom baryons are recognized to be 
one of the main research areas in particle physics. During the last few years quite interest- 
ing experimental observations have been obtained in heavy hadron spectroscopy. With the 
help of several refined measurements, the new states are discovered both in the charm and 
bottom sector (for recent experimental results see review [1]). These observations come 
from both the BaBar and BELLE, as well as from CDF and Dp Collaborations. LHC 
opens the possibilities for the discovery and detailed study of the new baryon states [2]. 
The considerable progress on the experimental side, has stimulated the theoretical investi- 
gation for understanding the dynamics of heavy flavor hadrons. Careful and comprehensive 
theoretical studies of the experimental results on heavy hadron spectroscopy and analysis 
of their weak and strong decays can provide essential knowledge on the quark structure of 
these hadrons. The strong coupling constants of the light pseudoscalar and vector mesons 
with heavy baryons are the main parameters for understanding the dynamics of heavy 
baryons. For this reason, reliable determination of the strong coupling constants of light 
pseudoscalar and vector mesons with heavy baryons within QCD receives special atten- 
tion. Unfortunately, at the hadronic scale, QCD becomes nonperturbative and it makes 
impossible to calculate these coupling constants starting from QCD Lagrangian. Therefore, 
for calculation of these couping constants some nonperturbative methods are needed. The 
QCD sum rule approach, which is based on the fundamental QCD Lagrangian, is one of 
the most attractive and applicable approaches [3]. In this work, we calculate the strong 
coupling constants of light vector mesons with heavy baryons within the light cone version 
of the QCD sum rules (LCSR) (for more about LCSR see [4]). The coupling constants of 
pseudoscalar mesons with heavy baryons is studied in detail in the same framework in [5]. 

The work is arranged in the following way. In section 2, the light cone sum rules for the 
coupling constants of light vector mesons with heavy baryons are obtained. In the following 
section, the numerical analysis of the obtained sum rules is performed and a comparison of 
our results with ones existing in literature is presented. 

2 Light cone QCD sum rules for the heavy baryon— 
light vector meson couplings 

In this section, we calculate the strong coupling constants of light-vector mesons with heavy 
baryons. Before starting to calculate these coupling constants few words about the SU (3) / 
classification of the heavy baryons are in order. Baryons with a heavy single quark belong 
to either antisymmetric antitriplet 3^ or symmetric sextet 6^ representations. Using the 
symmetry properties of the wave function, the spin of the light diquark is equal to zero for 
the antitriplet, while it is equal to one for the sextet. For this reason the total spin of the 
ground state baryons is 1/2 for 3f, but it can be both 3/2 and 1/2 for 6^. In the present 
work we restrict ourselves by considering only spin 1/2 heavy baryons. 

In order to obtain the strong coupling constants of light-vector mesons with heavy 
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baryons within the LCSR method, we consider the following correlation function: 



n fe) = % J d 4 xe ipx (V(q) |r {^(aO^O)}! 



(1) 



where the indices i and j get two values and describe the sextet-sextet (i = 1, j = 1), 
sextet-triplet (i — 1, j — 2) and triplet-triplet (i = 2, j = 2) transitions, respectively. 
In further respect, we introduce the definitions, nJ 11 ) = n^, nJ 12 ) = 11^ and nJ 22 ) = 
U^ 3 \ correspondingly. Note that, the V(q) in Eq. (1) denotes the light-vector mesons 
(p, u, K*, 0) with momentum q and rj is the interpolating current of the heavy baryon. 

The correlation function (1) can be calculated in two different kinematical regions, 
namely, in terms of hadrons (phenomenological side), as well as in deep euclidean region 
when p 2 — > — oo, in quark and gluon degrees of freedom using operator product expan- 
sion (OPE) (theoretical side). Equating both representations with the help of dispersion 
relations allows us to obtain the sum rules that will be used for calculation of the strong 
coupling constants of the light-vector mesons with heavy baryons. 

The phenomenological part of the correlation function can be obtained by saturating 
it with hadrons which carry the same quantum numbers as their interpolating currents. 
Separating contribution of the ground state baryons we get, 



(») 
%2 



B 2 (p)) (B 2 (p)V(q) | B 1 (p + q)) (B 1 {p + q) 



1 IB 1 



(p 2 — ml) [(p + q) 2 — ml] 



+ 



(2) 



where m x and m 2 are the masses of the initial and final baryons and dots represent contri- 
butions from higher states and continuum. The matrix elements appearing in Eq. (2) are 
defined as follows: 







Vb 



B(p)) = \ t u(p) , 



(B 2 (p)V(q)\ B 1 (p + q)) =u(p) 



m l + m 2 



u(p + q)e 11 , 



(3) 
(4) 



where Aj are the residues of the heavy baryons and are the momentum and vector 
polarization of the vector meson, fi and f 2 are the charge and magnetic form factors, 
respectively. 

Performing summation over spins of the baryons and using Eqs. (2)-(4), we get for the 
phenomenological part, 



jj(y') 



XiXj 



(p 2 - m 2 ) [{p + 



g)2 _ m 2j {(A + h)Hi + 2 fi( £ -p)^ + other structuresj . (5) 



The reason why we choose the structures p 1 ^ and (e-p)^> is that they show the best conver- 
gence. 

In calculating the theoretical part from the QCD side, explicit expressions of the inter- 
polating currents of heavy baryons are needed. Using the fact that the sextet (antitriplet) 
current should be symmetric (antisymmetric) with respect to the light quarks, the most 
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general form of the interpolating currents for the spin-1/2 sextet and antitriplet baryons 
can be written in the following form, 



(s) 1 



abc 



{ (qf CQ b )l 5 q c 2 - (Q aT Cqi) l5 q<i + [(qf C l5 Q b )q c 2 - (Q aT C l5 qi)ql\} , 



V ( q ] = ^f^^if Cq b ) l5 Q c + (qf CQ b ) l5 q c 2 + (Q aT Cq b 2 ) lb q{ 
+ p[2(qfC lb q b 2 )Q c + (qf C l5 Q b )q c 2 + (Q aT C im b 2 )q^) , 



(6) 



where the superscripts s and a refer to the sextet and antitriplet, respectively, and subscripts 
a, b, c are the color indices, C is the charge conjugation operator, j3 is an arbitrary parameter, 
and (3 = —1 case describes the Ioffe current. The light quark fields q\ and q 2 for the sextet 
and antitriplet are given in Table 1. 





91 


<?2 


v +(++) 
v 0(+) 

y-(0) 

w-(oy 

"6(c) 
~0(+)' 

^ b( (o) 

n 6(c) 

A o(+) 

"6(c) 

~0(+) 

"6(c) 


■u 

d 
s 

u 


u 
d 
d 
s 
s 
s 
d 
s 
s 



Table 1: The light quarks qi and q 2 for the sextet and antitriplet baryons 

Before calculating the correlation functions responsible for all sextet-sextet-vector mesons 
(SSV), sextet-antitriplet-vector mesons (SAV) and antitriplet-antitriplet-vector mesons 
(AAV) transitions from the QCD side, we find the relations among invariant functions for 
each of the above-mentioned classes (see also [5-10]). As a result, we find that the magnetic 
and electric couplings for each class of transitions are described in terms of only one invari- 
ant function. Of course, the invariant functions for each class of transitions are different 
from each other in the general case. It should also be noted that the relations among the 
invariant functions are all structure independent. 

After these remarks, we can proceed to establish relations among invariant functions 
involving the couplings of sextet-sextet transitions. Let us first consider S° — > S°p° tran- 
sition. The invariant function responsible for this transition can schematically be written 
as: 
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where we have introduced the formal notations, 

U { i\u,d,b) = (uu\E b E° b \0) , 
U^(u,d,b) = (5&|E°E°|o) , 

ni (1) (u,d,6) = (Jd|E°E°|o) . (8) 
The interpolating current of p° is formally written in the form, 

jf = Yl dpqqQ^Q , (9) 



u,d,b 



where we have set g p o bb = and g p o Uu = —g p o dd = l/y/2. Physically, each term on the right 
hand side of Eq. (7) describes emission of the p° meson from u, d and b quarks of the E° 
baryon, respectively. Since interpolating current of E° is symmetric under the exchange 
of u and d quarks, it leads to the result, n/^(u, d, b) = U^\d, u, b). Using the definitions 
given in Eq. (8) and taking the remark after Eq. (9) into consideration, we obtain, 



nyO,yO n O 1 

V2 



n [ i\u,d,b) -u?(d,u,b) . (io) 



Obviously, in the SU(2) f limit n E °^ E °"° = 0. 

The invariant function describing E^ — > S^p° (E^ — > E^p°) can be obtained from the 
invariant function for the E° — > E°p° transition with the help of the replacement, d — >■ u 
(u —> d), and using the fact that E° = -\/2T^ (y/2^). As a result, we obtain, 

4H { i\u,u,b) = -2{uu |S+E+| 0) , (11) 
4n^ 1} (d, d,b) = 2 (du | E~E- 1 0) . (12) 

Since E^ - -* contains two u(d) quarks, there are four possible ways for emitting p° from the 
u(d) quark. Hence, the related invarian functions are obtained as: 

n*t^tf>° = V2n?(u,u,b) , (13) 

n E ^ E ^° = Vm { i\d,d,b) . (14) 

The result for the invariant function responsible for the H^ " 1 — > E b ~^ p° transitions 
can easily be obtained from the result for the E° — > S°p° transition using the fact that 
= E°(d s) and Si" = E°( M s), i.e., 

U^'>° = -Lll?\u,s,b), (15) 
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U^^^' = -—Uf\d,s,b). (16) 

How can one find the relations among the invariant functions in the presence of the 
charged p* 1 mesons? In order to answer this question, we again consider the matrix element, 
(dd |E°E°| 0). This matrix element means that the d quarks from E° and E° form the final 
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dd state, while u and b quarks are the spectators. The matrix element (ud 0) 
corresponds to the case when d quark from E° and u quark from form the ud state, 
again the remaining u and b quarks being the spectators. This fact allows us to comment 
that these matrix elements are proportional to each other. A detailed calculation shows 
that, these matrix elements are related to each other through, 

I&^tr- = (ud |E+E°| 0) = -V2 (dd |E°E°| 0) 

= -V2U?(d,u,b) . (17) 

Making the replacement u <-y d, from Eq. (17) we get, 

n^^' + = Vzn?\u,d,b) . (is) 

In calculating the coupling constants of SSV, SAV and AAV, it is enough to consider 
the transitions E° — y T^V ', — y Sj,°V and S° — y E^V, respectively, since all other strong 
couplings can be achieved from these results with the help of corresponding replacements 
among the quarks. In order to obtain the correlation functions in terms of III, which 
describe transitions among sextet-sextet, sextet-antitriplet and antitriplet-antitriplet with 
other vector mesons, similar calculations can be done. These results are presented in 
Appendix A. 

We can now proceed to calculate the invariant functions (i = 1,2,3) responsible 
for E° — y S°p°, S 6 ° — y H 6 °p° and S° — y H°p°transitions from QCD side. These invariant 
functions can be calculated in deep Eucledian region — p 2 —y oo, — {p + q) 2 — > oo using 
the operator product expansion. The main nonperturbative ingredient in the calculations 
are the distribution amplitude (DA's) of the vector mesons. These distribution amplitudes 
appear in determination of the matrix elements of nonlocal operators (V(q) \q(x)Fq(0)\ 0) 
and (V(q) \q(x)G ^?(0)| 0), where T is any Dirac matrix. Up to twist twist-4 accuracy, the 
expressions for the distribution functions of vector mesons can be found in [11, 12]. 

In calculating the correlation functions responsible for afore-mentioned decays from 
QCD side, the expressions of light and heavy quark propagators are also needed. The light 
quark propagator in an external field is calculated in [13] whose expression is given as 



-ig s / du 



o 



(19) 



where — 0.577 is the Euler Constant, and A is the scale parameter and it is chosen as 
the factorization scale A = (0.5 -r 1) GeV (for more detail see [14]). 
The propagator for the heavy quark is [15]: 



m 2 Q KrirriQ^ 2 ) . m 2 Q £ 

^ ) = 4VWZP *i^^(m V^) 

" Ws I T^ze-^ f du\ ^ + mQ — ' U 



o 



G^(ux)a^ + — u x,G^ lv 



2{m 2 Q — k 2 ) 2 ^ u m, 



, (20) 
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where Ki are the modified Bessel function of the second kind. 

Using the expressions of light and heavy quark propagators, as well as definitions of the 
DA's for the vector mesons, and after lengthy calculations one can obtain the correlation 
function from QCD part. 

As has already been noted, the relations among the correlation functions for the con- 
sidered transitions are structure independent, but their explicit expressions are structure 
dependent. For this reason we introduce new indices a in the correlation function, where 
a = 1 stands for the choice of the structure (s-p)]t> and a = 2 for the structure j^. 

Equating the coefficients of the structures, (ep)^> and fyfi for the hadronic and QCD sides 
and performing Borel transformations over the variables p 2 and (p + q) 2 , which suppress 
the contributions of the continuum and higher states, we finally get the sum rules for the 
strong coupling constants of light vector mesons with sextet and antitriplet heavy baryons 
as: 

(i)2 (i)2 2 
1 m l , m 2 , m V 

f(») — P M? ~mT M( + M% TTW (Q-\\ 

where M 2 and Mf are the Borel parameters corresponding to the initial and final heavy 
baryons, respectively. In the problem under consideration the masses of the initial and final 
heavy baryons are very close to each other, hence we can take Mf = Mf = 2M 2 . Note 
that the residues of the sextet and antitriplet heavy baryons are calculated in [16]. This 
leads to the result that for numerical calculation of the sum rules, the DA's are needed to 
be evaluated only at u = M ^ M s = \- 



3 Numerical analysis 

Having already obtained the sum rules for the strong coupling constants of light vector 
mesons with heavy baryons, we can now proceed evaluating them numerically. The essential 
ingredients of the LCSR are the DA's of the light vector mesons. Explicit expressions of 
DA's for the vector mesons and the values of the parameters entering to the expressions of 
DA's are given in [11, 12]. The residues of the heavy baryons are calculated in [16]. 

In addition to the DA's and other input parameters, the sum rules for SSV, SAV and 
AAV transitions contain three auxiliary parameters : the continuum threshold s , Borel 
parameter M 2 and the parameter /3 entering the expressions of the interpolating current. 
For this reason we try to find such regions of these parameters where strong coupling 
constants are practically independent of them. 

In finding the working region of M 2 , we require that the continuum and higher state 
contributions should be less than half of the dispersion integral, and additionally the contri- 
bution of the higher terms with the power 1/M 2 be 25% less than the total result. These two 
restrictions lead to the result that the "working region" of M 2 is 15 GeV 2 < M 2 < 30 GeV 2 
(for the bottom baryons) and 4 GeV 2 < M 2 < 12 GeV 2 (for the charmed ones). The con- 
tinuum threshold is varied in the region (m B + 0.5) 2 GeV 2 < s < (m B + 0.7) 2 GeV 2 . 
To be more illustrative about how the numerical analysis is performed, as an example, we 
consider the S b ° — > S b °p° transition. In Figs. (1) and (2) we present the dependence of fi 
and j'2 on M 2 for the above-mentioned transition, at five different values of (3, and at a 
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Jl 


5.0±1.7 
10.0±3.4 


2.6±0.9 
7.0±2.4 


'4- '4- 


4.0±1.3 
11.0±4.0 


2.5±0.8 
23.0±8.0 



Table 2: The values of the strong coupling constants f\ for the transitions among the 
sextet-sextet heavy baryons with vector mesons. 

fixed value of sq. From these figures, one can conclude that the strong coupling constants, 
fi and f 2 are practically independent of M 2 when it varies in its own 'working region" . 

In Figs. (3) and (4) we present the dependence of f\ and f'2 on cos#, at three fixed 
values of sq and M 2 , where tan# = (3. We observe from these figures that when cos 6? is 

varied in the region —0.5 < cos# < 0.3, the results are insensitive to the variation of fi. 

s' -»s'°p s' ->s'°p 
From these figures, we obtain, f~ b b H = 2.2 ± 0.7 and fp b = 30 ± 10. Similar 

analysis for the other couplings of vector mesons are carried out and the results for j\ and 

f 2 are presented in Tables (2)-(7), respectively. The errors presented in these Tables are 

due to the variation of the auxiliary parameters, as well as uncertainties in the values of 

the input parameters. 

The results on the coupling constants of the light vector mesons with heavy baryons 
presented in Tables (2)-(7) lead to the following conclusions: 

• For the coupling constant, f\: The predictions for the strong coupling constant f\ 
which are obtained using the most general and Ioffe currents (/3 = —1) disagree 
considerably from each other, especially for the channels, S 6 ° — > S c + — > E+lu, 

S+ S>, ~' b ° -+ E~u, S- Kip- and £° -> A+p". 

• For the /2 channel: predictions of the general and Ioffe currents for the strong coupling 
constants of SSV and AAV show considerable differences. These discrepancies can be 
attributed to the fact that, (5 = — 1 lies outside the stability region of f3 and it makes 
the predictions less reliable. 

Our final remark in this section is that the coupling constants of the S° — > H°p° , 
E° c — > E° c uo and E° c — > E° c 4> transitions are also calculated in the HQET [17] within 
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Table 3: The values of the strong coupling constants f\ for the transitions among the 
sextet-antitriplet heavy baryons with vector mesons. 
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4.9+1.6 
4.6+1.5 


4.3+1.4 
4.1+1.4 



Table 4: The values of the strong coupling constants f\ for the transitions among the 
antitriplet-antitriplet heavy baryons with vector mesons. 

the same framework as our work. Our predictions on the coupling constants of these 
transitions and those given in [17] are in a good agreement with each other. 
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Table 5: The values of the strong coupling constants fa for the transitions among the 
sextet-sextet heavy baryons with vector mesons. 

Appendix A : 

Here in this appendix we present the expressions of the correlation functions in terms 
of invariant function II ^ involving p, K, omega and mesons. 

• Correlation functions responsible for the sextet-sextet transitions. 
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Table 6: The values of the strong coupling constants fi for the transitions among the 
sextet-antitriplet heavy baryons with vector mesons. 
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Table 7: The values of the strong coupling constants fi for the transitions among the 
antitriplet-antitriplet heavy baryons with vector mesons. 



U n b -^ B b 

n 6 f> 



E, - »E K*° 



]JH ^b 



E+^E+ W 



n^6 



E-^-srw 



U?(d,u,b) , 
Tl?\d,u,b) , 

-^[nf ) («,d,6) + ni 1) (d,«,6) 

Vm ( i\d,d,b) , 
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rp^ = i=n« ( U ,s,b), 
n<"-^" w = -^n^d,^) , 

u s' b -^s'-4> = u^\s,d,b) , 

11^^ = 2U^\s,s, b) . (A.l) 
• Correlation functions responsible for the sextet-antitriplet transitions. 



nS ;-^ 6 v> = __L n ( 1) (rf, s ,6) ! 

n s b -^A 6P - = y/2n0-\ u , d,b) , 
n s;-^s-p- = n W(d, s ,6) , 

n H^H- P + = n( 1} ( M , s ,6) , 
n=8^ = -n< 1 )(d,«,6), 

n E 6-^6* ,0 = -v^n5 1) (d,d,6), 
n ^ s -/?- = y2n«( s , s ,6), 

n^^ = -n«(d, w ,fo), 

Tl^^ K *° = -V2Tl£\d,d,b) , 

n Q ;^ K *° = V2n£\ s ,s,b) , 
u<°^ K "° = -n£\d,u,b) , 
n^A 6 x*+ = -V2ii£\u,u,b) , 

T&^k** = -I$\ u ,d,b), 
n n fc --sgA:- = V2n«( S , S ,6) , 

n s l-^-=n«( M , s ,6), 
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U?(u,d,b)-U^>(d,u,b) 



1 

-n$-\ s ,u,b) 

-I$\s,d,b) 



• Correlation functions responsible for the antitiplet-antitriplet transitions. 



(A.2) 



n=6 



-^=uf\u,s,b), 



1 



-—U^(d,s,b) , 



j-jA 6 ^A 6/ 9° }_ 



1 

7^ 



n? ) (d,«,6)-n^(«,d,6) 



a), 



n ag-.sg w 



n«( M , s ,6) , 
n«( M)M) 6) , 
nW(«,«,6) , 

i 
1 

72 

1 



6 ~ — 



v/2 



(A.3) 



The expressions for the charmed baryons can easily be obtained by making the replace- 
ment b — > c and adding to charge of each baryon a positive unit charge. 
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Figure 1: The dependence of the f{ b ~ 6 on M 2 at fixed values of the s and 
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Figure 2: The dependence of the ^ 6 on M 2 at fixed values of the s and /3. 
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